Different escape modes in two-photon double ionization of helium 
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The quadrupole channel of two-photon double ionization of He exhibits two distinctly different 
modes of correlated motion of the photoelectron pair. The mode associated with the center-of- 
mass motion favors a large total momentum which is maximazed at a parallel emission. However, 
the mode associated with the relative motion favors an antiparallel emission. This difference is 
manifested in a profoundly different width of the angular correlation functions corresponding to the 
center-of-mass and relative motion modes. 
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The process of correlated motion of multiple ionization 
fragments has been at the forefront of atomic collision 
physics during the past decade. Recent progress in exper- 
imental techniques made it possible to detect simultane- 
ously a large number of charged reaction fragments with 
fully determined kinematics [l| . The long range Coulomb 
interaction between these fragments makes a full theo- 
retical description of such a process a highly challenging 
task. In the meantime, the simplest multiple fragmenta- 
tion reaction, the single-photon double ionization (SPDI) 
of helium is now well understood with accurate theoret- 
ical predictions being confirmed experimentally under a 
wide range of kinematical conditions iii. All the 
information about the correlated motion of the photo- 
electrons is described in SPDI by a pair of symmetrized 
amplitudes f ± {9\2, E\,E 2 ) which depend on the relative 
interelectron angle and energy 0, [g. If the photon is 
described by the linear polarization vector e, the dipole 
matrix element of SPDI can be written simply as 



D oc /+ (fer + k 2 ) ■ e + f~ (fei - fca) • e 



(1) 



where ki = fej/fcj, i = 1,2 are the unit vectors directed 
along the photoelectron momenta ki. Generalization of 
Eq. (D to arbitrary polarization of light is straightfor- 
ward pfl. Under the equal energy sharing condition, the 
anti-symmetric amplitude vanishes f~(E\ = E 2 ) = 
and all the information about the SPDI process is con- 
tained in one symmetric amplitude / + . Following pre- 
dictions of the Wannier-type theories H, 0, the SDPI 
amplitude can be written using the Gaussian ansatz: 



|/+| 2 oc exp 
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where the width parameter A#i2 indicates the strength 
of angular correlation in the two-electron continuum. Al- 
though the analytical theories 0, 0] validate Eq. (01 only 
near the double ionization threshold, numerical models 



[Tof and direct measurements 0, 113 su pport its validity 
in a far wider photon energy range [lCj. 

Two-photon double ionization (TPDI) of He is a much 
more complex fragmentation process with two compet- 
ing decay channels into the S and D two-electron con- 
tinua. In a most general case, complete separation of 
the dynamic and kinematical variables in TPDI requires 
introduction of five symmetrized amplitudes [Iflj or four 
unsymmetrized amplitudes |14|. 

In this Letter, we demonstrate that complexity of the 
TPDI leads to a new phenomenon of two distinct corre- 
lated escape modes: one is associated with the center- 
of-mass motion of the photoelectron pair whereas an- 
other is related to their relative motion. The center-of- 
mass motion favors a large total momentum of the pair 
p = k\ + fc2 which is gained at a parallel escape. In this 
configuration, the intcr-clcctron repulsion is strongest 
and the angular correlation factor has a relatively nar- 
row width. On the contrary, the relative motion favors 
a large relative momentum k — k\ — ki which is max- 
imized at an antiparallel escape. In this configuration, 
the inter-electron repulsion is weak. This results in a 
profoundly large angular correlation width as compared 
with the center-of-mass motion. 

For simplicity, we consider equal-energy-sharing kine- 
matics Ei = Ei and restrict ourselves with only the 
dominant quadrupole TPDI amplitude. By performing a 
derivation similar to Ref. [13| , the quadrupole amplitude 
can be parametrized as 

Q oc g k {k ® k} 2 +g P {p®p}2 (3) 
+£/o|[fe x p] <8> [k x p]| • {e ® e} 2 

Here we introduced the unit vectors p = p/p and k = 
k/k. The symmetrized amplitudes in Eq. © are defined 

as 



9k,p 



(4) 
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±2Pl' +1 {x)]Qu +2 {ki,k 2 ) T -CiP;'(x)Q ll (k 1 ,k 2 ), 
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g Q (x) = ^2C l Pl'(x)Q ll (k u k2). 
1=1 

Here the upper and lower sets of signs refers to gk and 
g p , respectively. These amplitudes depend on the in- 
terelectron angle x = cos 9 12 and are expressed via the 
quadrupole radial matrix elements Qw(ki, fe) and the 
Legendre polynomial derivatives P('(x). Normalization 
coefficients Ci and Di are given in Ref. Using sim- 
ilar notations, the dipole amplitude of SPDI Q under 
the equal energy sharing condition is parametrized as 
D oc f p (p-e), /+ = f p , f~ = 0. We note that D is linear 
with respect to p and does not contain k under the equal 
energy condition. In contrast, Q is quadratic with respect 
to k and p and contains both vectors even when E\ = E2 ■ 
The amplitude gk which enters Eq. © with the tenso- 
rial product {k <gi fc} 2 can be associated with the relative 
motion of the photoelectron pair described by the vector 
k. Similarly, the amplitudes g p can be associated with 
the center-of-mass motion and the amplitude go which is 
entering Eq. © with the vector product k x p can be 
associated with the mixed motion mode. 
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FIG. 1: The TPDI amplitudes in the quadrupole channel 
gk (top) and g p (bottom) fitted with the Gaussian ansatz 
The excess energy of 4 eV is shared equally between the 
photoelectrons E\ — E2 — 2 eV. The arrows indicate the 
Gaussian width parameter A#i2. 

To calculate the radial matrix elements Qw(ki,k2), 
we employed here the same dynamical model as was 
outlined in our previous work |13| . In this model, the 
electron-photon interaction was treated in the lowest- 
order perturbation theory using the closure approxima- 
tion whereas the electron-electron interaction was in- 
cluded in full using the convergent close-coupling (CCC) 



method. The model proved to be capable of describing 
the angular correlation pattern in the two-electron con- 
tinuum in good agreement with non-perturbative, with 
resp ect to the electromagnetic interaction, calculations 

We calculated amplitudes Q in a range of excess en- 
ergies Ei + E2 from 1 eV to 20 eV above the double 
ionization threshold. We employed a fairly large CCC 
basis set composed typically of 25 — I box-space target 
states 01 with < I < 6. Convergence of the calculation 
with respect to the basis size was thoroughly tested. 

In the whole excess energy range, the amplitude go 
was found insignificant as compared with gk and g p . The 
latter amplitudes were fitted with the Gaussian ansatz 
. A typical quality of the fit can be judged from Figure 
n where the amplitudes gk and g p are exhibited for E\ — 
E2 = 2 eV. The corresponding width parameters A#i2 
are plotted in Figure |2 as a function of energy, along 
with the width parameter of the dipole amplitude f p . 
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FIG. 2: The Gaussian width parameters A612 of the am- 
plitudes <7fc (red filled circles), g p (blue open circles) and f p 
(light blue squares) as functions of the excess energy E\ +£2- 
Extraction of the width parameters is illustrated in Figure Q 
for Ei = E 2 = 2 eV. 

We observe a much larger Gaussian width of the rela- 
tive motion amplitude gk as compared with the center-of- 
mass motion amplitudes g p and f p , the latter two having 
a very similar width. We interpret this stark difference 
in terms of the strength of the electron-electron repul- 
sion. This strength is much larger in the center-of-mass 
motion which favors large p and hence parallel emission, 
as opposed to the relative motion which favors large k 
associated with anti-parallel emission. We note that the 
mixed mode amplitude is typically an order of magnitude 
smaller than the pure mode amplitudes gk and g p thus 
supporting our notion of distinct escape modes. 

In the case of co-planar geometry when all three vec- 
tors fci , k 2 and e belong to the same plane, we can choose 
the coordinate frame as z || ki and x \\ [ki x fea]. In this 
case, Eq. {Q and Eq. © are simplified to 

D cx f p (cos 81 + cos 82 ) (5) 
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(cos 6*i — cos# 2 ) 2 — -(1 - cos 9i 2 ) 



(COS 0\ + COS 62 ) --(1 + COS #12) 



: 9o (cos 



1) 



(6) 



From these equations, it is seen that the amplitudes f p 
and g p related to the center-of-mass mode, and the ampli- 
tude gk related to the relative motion mode, contribute 
quite differently to the corresponding matrix elements. 
All three amplitudes peak strongly near #12 = 180°. 
However, the kinematic factors corresponding to f p and 
g p have nodes at this angle whereas the kinematic factor 
accompanying gk has a peak. As a result, the term pro- 
portional to gk dominates strongly the quadrupole ampli- 
tude. This dominance is illustrated in Figure|31where the 
triply-differential cross-section da / dEidflid^ oc \Q\ 2 of 
the TPDI of He at E x = E 2 = 2 eV and the coplanar 
kinematics is plotted as a function of the variable escape 
angle 2 - Various panels from top to bottom correspond 
to fixed escape angles X = 0°, 30°, 60° and 90°. Two 
calculations are displayed in the figure: one with the full 
quadrupole amplitude Eq. an d another with only the 
gk contribution. The dominance of the gk term is evi- 
dent at all fixed angles # 2 . It is particularly strong at a 
fix escape angle Q\ = when nearly all the contribution 
to the quadrupole amplitude comes from the gk term. 

We note that in our previous paper [l3| we introduced 
a slightly different set of amplitudes in the quadrupole 
channel: 



Q cx -^+[P2(cos0i)+P 2 (cos0 2 ) 



(7) 



3 COS 0\ COS 02 — COS 0\2 



3o(cos 2 6»i2 - 1) 



where 



9s 



9 P + 9k 
2 (&> - 9k) 



(8) 



Although the set of amplitudes © gives an identi- 
cal quadrupole amplitude, it does not provide such a 
clear separation of the center-of-mass and relative mo- 
tion modes. Our earlier attempt in Ref. [T^ | to apply the 
Gaussian ansatz J2J) showed no such a clear systematic 
behavior of the angular correlation width with respect of 
the excess energy as exhibited in Figure 
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FIG. 3: Triply-differential cross-section da / dEid£lid£l2 oc 
\Q\ 2 of the TPDI of He at Si = E 2 = 2 eV in the coplanar 
kinematics is plotted as a function of the variable escape angle 
62- Various panels from top to bottom correspond to fixed 
escape angles 6»i =0°, 30°, 60° and 90°. The red solid line 
corresponds to the full quadrupole amplitude © whereas the 
blue dotted line exhibits the contribution of the gt term only. 
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In conclusion, we demonstrated the presense of the 
two distinct photoelectron escape modes following two- 
photon double ionization of He. One of the modes corre- 
sponds to the center-of-mass motion of the photoelectron 
pair. It enhances the total momentum of the pair and 
therefore favors the parallel emission. The inter-electron 
repulsion is strong in this mode and the angular corre- 
lation width is relatively small. The other, relative mo- 
tion, mode enhances the relative momentum of the pair 
and therefore favors the antiparallel emission. The inter- 
electron repulsion is much weaker in this mode and the 
angular correlation width is march larger. Both modes 
are fully symmetric and present under the equal energy 
sharing condition. In contrast, the single-photon double 
ionization has only one fully symmetric mode which is 



associated with the center-of-mass motion. This mode, 
in terms of the angular correlation width, is very sim- 
ilar to the center-of-mass motion mode in two-photon 
double ionization. The presense of two modes is sim- 
ply a reflection of the quadratic tensorial structure of 
the quadrupole photoionization amplitude as compared 
to the linear structure of the dipole photoionization am- 
plitude. 
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